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摘 要 : 研究 典型 株 型 沙 生 治 从 周围 的 流 场 分 布 , 旨 在 为 干旱 、 半 干旱 地 区 合理 选择 不 同 株 型 的 防风 沙 植被 提供 
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J FLUENT 软件 对 3 类 典型 株 型 ( 坛 形 、 梭 形 \ 帅 形 ) 灌 从 周围 的 流 场 进行 数值 模拟 ,分 析 不 同 植 


株 形态 对 风沙 流 的 影响 ,并 加 以 风 洞 试验 验证 ,结果 表明 :(1) 3 类 株 型 周围 流 场 可 分 为 5 个 区 , 且 株 后 均 存 在 3 个 涡 
流 。 受 涡流 强度 的 影响 ,在 积 沙 初始 阶段 , 梭 形 . 帅 形 植株 主要 在 株 后 6~7 HH 处 积 沙 ,而 坛 形 植株 在 3 HH 附近 积 沙 。 
(2) 受 植 株 最 大 侧 影 面积 高 度 层 的 影响 ,3 类 株 型 灌 丛 株 后 1 五 处 的 风速 极 小 值 依次 出 现在 0.3 m、0.4 m、0.8 m 高 度 
处 ,最 优 防护 高 度 依次 为 0.2~0.4 m、0.3~0.6 m、0.8~1 m。3 类 株 型 株 后 的 空气 动力 学 粗糙 度 逐 渐 减 小 , 且 坛 形 的 粗 
糙 度 明显 高 于 其 他 株 型 。(3) 3 类 株 型 在 -2~10 HH 范围 内 均 可 有 效 降低 风速 , 株 后 近 地 表 区 防风 效益 表现 为 坛 形 > 
梭 形 > 帅 形 ,而 中 高 空 区 防风 效益 均 随 株 距 增 加 而 减 小 ,(4) 在 ?=10s 时 ,3 类 植株 周围 总 积 沙 长 度 分 别 为 8.5 H6 
H、4.5 H, 梭 梭 、 沙 拐 京 分 别 在 距 入 口 5~5.5 m、4.5~6 m 处 存在 不 同 程度 的 风 刨 现象。 对比 其 他 植株 , 白 刺 (Nitraria 
sphaerocarpa ) 具 有 较 好 的 阻 沙 效果 ,在 防风 固沙 工程 建设 中 建议 将 其 与 梭 梭 (Haloxylon ammodendron ) WPR ( Cal- 


ligonum mongolicum ) 结 合 , 赋 能 发 挥 白 刺 的 阻 沙 作 用 ,又 可 利用 梭 梭 、 沙 拐 束 较 好 的 中 高 空 防风 效果 。 


关键 词 : 沙 生 灌 从 ; 植物 株 型 ， 防风 固沙 ; 数值 模拟 


植被 对 地 表 起 到 保护 作用 ,增加 植被 覆盖 率 是 
防风 固沙 的 主要 手段 "。 然 而 ,在 干旱 、 半 干旱 地 
区 ,由 于 水 分 条 件 的 限制 ,植被 主要 以 低 盖 度 沙 早 
生 灌木 为 主 , 许 多 灌木 以 单 株 或 单 丛 的 形式 存在 ， 
其 与 风蚀 关系 的 研究 可 简化 为 单 植株 的 防风 周 沙 
研究 2”。 但 不 同类 型 植株 对 周围 气流 场 的 影响 存在 
差异 ”。 因 此 ,研究 不 同 株 型 单 株 植被 周围 的 风速 
与 输 沙特 征 , 对 干旱 、 半 干旱 区 合理 选择 风蚀 防治 
中 的 植被 类 型 具有 重要 意义 ,同时 也 对 建设 不 同形 
态 的 稀 疏 灌 从 防护 林 工 程 具有 一 定 的 参考 意义 。 

株 型 会 显著 影响 植株 的 防风 固沙 效益 ,目前 有 
关 植 株 形态 与 风蚀 关系 的 研究 主要 存在 两 种 方法 ， 
一 类 是 采用 传统 的 风 洞 试验 或 野外 试验 探究 植株 
形态 对 风蚀 影响 的 差异 ,如 亢 力 强 等 “通过 风 洞 试 
验 测 量 了 细 长 状 植株 和 上 大 下 小 形状 植株 对 空气 
动力 学 粗糙 度 的 影响 ;Abbas 等 ”对 两 种 株 型 植株 进 
行 研 究 , 指 出 在 不 同 风 速 条 件 下 , 女 上 页 (Ligusirum lu- 
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cidum) 对 强风 的 抵抗 力 大 于 秋英 (Cosmos bipinna- 
tus)。 此 类 研究 主要 探究 了 单 株 植物 对 周围 风 场 特 
征 的 影响 ,对 于 沙 粒 的 运动 及 植株 周围 输 沙 特征 关 
注 较 少 。 男 一 类 是 利用 数值 模拟 探究 单一 植被 要 
素 周 围 的 流 场 变化 , 相 比 野外 试验 或 风 洞 试验 , 数 
值 模拟 在 数据 获取 及 工作 量 等 方面 具有 明显 优 
势 。 李 正 农 等 "将 豆瓣 黄杨 LBuxrs sinica) 简 化 为 
3D 对 称 模型 ,有 效 模拟 了 树木 对 于 风 场 的 影响 ;Liu 
等 "通过 CFD 软件 研究 了 不 同形 状 单 株 植物 周围 的 
流 场 ,指出 底部 大 而 项 部 小 的 植物 具有 更 好 的 防风 
效率 。 以 上 研究 增强 了 人 们 对 植株 形态 影响 单 株 
植被 防风 效果 的 认识 ,但 在 现 有 单一 植被 要 素数 值 
模拟 中 ,更 侧重 探究 乔木 等 高 大 植株 周围 的 风 场 变 
化 ,对 于 治 从 及 不 同 株 型 灌 从 周围 风速 与 输 沙特 征 
的 研究 较 少 。 

沙 生 灌 从 在 沙 质 苇 漠 生态 系统 中 占有 重要 地 
位 ,维持 了 干旱 、 半 干旱 地 区 生态 环境 的 稳定 ”。 沙 
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生 灌 从 根系 发 达 , 生 合力 硕 强 ,在 生长 季 枝 叶 密 度 
较 大 ,区 域 植 被 盖 度 相应 增加 ,可 有 效 降 低 风 动量 ， 
从 而 起 到 防风 阻 沙 作 用 。 植 被 形态 和 结构 是 控制 
风蚀 的 优先 参数 ,不 同 株 型 灌 从 的 防风 阻 沙 效 
益 存 在 显著 差异 。 本 文 基于 数值 模拟 ,以 3 种 典型 
株 型 的 荒漠 汐 丛 为 研究 对 象 ,通过 分 析 不 同类 型 灌 
从 前 后 风速 .空气 动力 学 粗糙 度 ` 沙 粒 体积 分 数 等 
变化 规律 ,探讨 不 同 株 型 植被 的 防风 阻 沙 效 益 ,对 
比 得 出 优势 灌 从 株 型 ,并 利用 已 有 风 洞 试验 进行 验 
证 , 旨 在 揭示 不 同 植株 形态 对 风沙 流 活 动 的 影响 存 
在 差异 ,促进 干旱 、 半 干旱 地 区 防风 固沙 工程 设计 
中 对 灌 从 株 型 的 重视 ,为 防风 固沙 工程 及 生态 恢复 
工程 中 合理 选择 不 同 株 型 的 防风 沙 植被 提供 参考 。 


1 数据 与 方法 


1.1 研究 对 象 

本 文 研究 区 位 于 马兰 布 和 沙漠 东北 部 ,区 域内 
生长 着 不 同 株 型 的 荒漠 沙 生 灌 从 ,以 坛 形 RIE ir 
形 居 多 ,如 白 刺 (CNitraria sphaerocarpa ) Ki 4% ( Cara- 
gana Korsjinsjz)、 霸 王 (Sarcozygium Bunge) 、 杨 柴 


(Hedysarum mongolicum ) , > AŒ (Calligonum mon- 
golicum ) RIR (Haloxylon ammodendron ) 、 花 棒 ( He- 
dysarum scoparium) SE ,对比 选取 白 刺 RRAK 
沙 拐 束 作为 3 种 典型 株 型 参照 。 模 拟 所 需 植被 特征 
参数 源 于 已 有 野外 调查 试验 ” ,调查 于 灌 从 生长 季 
进行 ,在 林场 灌 从 均匀 分 布地 带 设 置 3 个 100 mx 
100 m 的 样 方 ,抽取 20 株 高 度 在 1 m 左 右 的 同龄 灌 
木 ,记录 高 度 、 冠 幅 等 植被 特征 参数 ,计算 平均 值 作 
为 确定 3 类 株 型 模型 的 参考 依据 ,结果 如 表 1 所 


示 。 本 文 主要 探究 二 维 尺 度 下 植株 周围 的 流 场 特 
征 ,为 便于 植株 建 模 , 对 植株 高 度 、. 冠 幅 进行 均值 处 
理 及 整数 化 处 理 。 经 过 试 算 与 模拟 ,确定 3 类 株 型 
模型 植被 参数 如 表 2 所 示 。 

由 表 1 可 知 , 白 刺 的 最 大 侧 影 面积 出 现在 0.10~ 
0.30 m 人 处 ,呈现 底部 大 、 顶 部 较 罕 的 坛 形 , 枝 下 高 度 
较 低 ,整体 贴近 地 面 ,重心 偏 低 ; 梭 权 整 体 呈 现 中 间 
HAKE IRIE ,形态 较为 高 大 ,重心 在 植株 中 间 ; 
沙 拐 囊 的 枝 下 高 度 最 大 ,最 大 侧 影 面积 出 现在 植株 
上 部 ,呈现 帅 形 ,重心 偏 高 。 
1.2 数值 模拟 方法 
1.2.1 植株 建 模 ”由 于 沙 粒 在 风沙 流 运 动 中 所 受 的 
RJ .重力 基本 在 同一 平面 内 , 故 本 文采 用 二 维 简 
化 模型 。 根 据 表 2 及 前 人 建 模 经 验 ,采用 GAM- 
BIT 软件 ,充分 利用 灌 丛 形态 的 对 称 性 ,建立 了 白 
刺 RR 沙 拐 吏 的 孔隙 几何 模型 ,依次 为 坛 形 AR 
FE . 灵 形 。 简 化 植株 模型 是 植株 形态 的 近似 和 概 
化 ,通过 较 小 的 植株 孔隙 度 以 及 各 类 实测 植被 特征 
参数 来 表现 生长 季 灌 从 的 植株 形态 ,因此 ,模拟 结 
果 具 有 一 定 的 可 靠 性 。 通 过 多 次 试 算 和 模拟 ,流域 
长 度 设 置 为 30 m, 计 算 高 度 为 10 m, 沙 床 厚度 为 
0.05 m; 植 株 简 化 模型 距 入 口 5 m, 计 算 域 简化 图 如 
图 1 所 示 。 

模型 形状 较为 规则 , 故 网 格 划 分 类 型 采用 Quad 
形式 (四 边 形 网 格 ), 网 格 划 分 方法 采用 结构 化 网 格 
划分 。 本 文 主要 讨论 沙 粒 的 流 场 特征 以 及 植株 前 
后 沙 粒 和 沙 床 表面 的 相互 作用 ,因此 ,对 地 面 及 植 
株 附 近 的 网 格 进行 加 密 ,经 过 多 次 模拟 和 调整 , 确 
定 网 格 尺寸 增长 率 为 1.04, 计 算 域 网 格 总 数 为 
44647 个 ,最 小 正 交 质量 远 远 大 于 0.99 ,最 大 正 交 算 


表 1 已 有 野外 调查 结果 


Tab. 1 Field survey results available 


植被 名 称 高 度 /m 冠 幅 /m 枝 下 高 度 /m 最 大 侧 影 面积 高 度 层 /m ”植物 侧 影 孔隙 度 /% ”植株 形态 
白 刺 0.834 0.106~0.828 0.03 0.10~0.30 34.54 坛 形 
TRIR 1.100 0.906~0.718 0.08 0.30~0.60 54.93 梭 形 
WPR 1.040 0.940~0.578 0.12 0.60~1.00 61.08 mI 
表 2 植株 模型 参数 
Tab.2 Model plant parameters 
植被 名 称 高 度 /m 迎风 面 冠 幅 /m 枝 下 高 度 /m 最 大 侧 影 面积 高 度 层 /m ”植物 侧 影 孔隙 度 /% 植株 形态 
白 刺 0.86 0.96 0.03 0.10~0.30 41.68 坛 形 
梭 梭 1.08 0.88 0.08 0.30~0.60 57.34 梭 形 
WPR 1.03 0.84 0.12 0.60~1.00 60.55 m 
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5m 固体 无 滑 移 边界 
(b) 梭 梭 WER 30m 


固体 无 滑 移 边界 
5m 固体 无 滑 移 边界 


人 对 称 边 界 30m 


0.84 m 
dao 着 0.12 m 
口 Ia 0.12m 
to 出 O 
: G 0.12 m 
固体 无 滑 移 边界 


固体 无 滑 移 边界 


图 1 流 场 简化 图 
Fig. 1 Model diagram 


斜率 远 远 小 于 0.001 ,网 格 质量 良好 ,满足 计算 要 求 。 
1.2.2 参数 设置 ”结合 前 人 "模拟 参数 设置 ,本 文 
将 风沙 流 视 为 不 可 压缩 流体 。 模 型 左 侧 入 口 处 采 
用 速度 边界 条 件 , 右 侧 采用 压力 出 口 边界 条 件 , 上 
壁面 采用 对 称 边界 ,下 壁面 及 植株 模型 采用 固体 无 
滑 移 边界 。 风 沙 流 粒 径 一 般 为 0.075~0.25 mm , 沙 粒 
粒 径 设置 为 4=0.1 mm, 沙 粒 密度 p=2650 kg* mi, 初 
台 沙 粒 体积 分 数 为 0.02% , 沙 床 初始 沙 床 堆积 率 a= 
0.625 , Hi u.=0.047 Pas, 空气 密度 p=1.225 kgm”, 
空气 运动 黏度 4=1.7894x10* Paso A OWA EEN 
典型 风速 廓 线 流 如 下 式 所 示 : 
v())= rar (1) 
式 中 : w AEE BARE Cms); yo WAREKE (m); k 
为 汉 卡 门 系 数 , 取 值 为 0.4;y 为 高 度 (m); vy) Ay Xb 
的 风速 值 (ms ')。 
求解 模型 采用 欧 拉 双流 体 模型 ,附加 汕 流 模 
型 。 由 于 需要 考虑 流 场 中 速度 变化 .求解 精度 问 
题 , 故 采用 非 定常 瞬 态 求解 方法 来 模拟 , 流 场 求解 
采用 Phase Coupled SIMPLE 算法 ,空间 离散 格式 采 
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(a) 白 刺 对 称 边界 ”30m 用 二 阶 迎 风格 式 , 时 间 步 长 取 0.001 s。 


1.2.3 基本 方程 ”本文 模拟 气流 为 不 可 压缩 流 , 主 
要 包含 的 控制 方程 有 连续 方程 .动量 方程 和 潮流 模 
型 方程 ,根据 模拟 需要 , 消 流 模型 选择 适合 高 雷诺 
数 的 标准 太 s 模 型 。 


连续 性 方程 为 : 
op a(pu,) a(p u,) _ 
ot + ðx i oy =o (2) 


AP: wu, u, 分 别 为 两 个 方向 的 速度 矢量 ; p 为 密 
度 (kg'm”); 为 时 间 (s)。 


动量 方程 为 : 
A(pu,) —\_ Op OT 。 aT, 
0 atat ay P 
a(pu,) E dp Ty T, 
FP +V(pu,@)=- 50+ a ae +p'g (4) 


式 中 : p 为 流体 微 元 体 上 的 压强 (Pa) ;五 为 速度 矢 
量 ; Ri Tosa WE pe Op 都 是 因 分 子 黏 性 作用 而 产生 
的 作用 在 微 元 体 表面 上 的 理性 应 力 7 的 分 量 (Pa); 
g 为 重力 加 速度 (ms?)。 

标准 太 s 模 型 的 庙 动 能 以 及 耗 散 率 方程 为 : 


M, 
| A 


k 


2 


de_ ð H, \ðe E e 
P dt ~ aC ae] = Cs AG + C,,6;) Csp 


(6) 
ee E pc, k Sie 
SUE: MAER ŽI u, = , C, NAME, C, 
RIR FA FKP RE BB EES IIE ; C, 表示 由 于 
FF FT 2 Me] | Go ES) ing 2 FE 5 i 为 时 间 (s); C,, =1.44; 
C,, =1.92; C,, =0.09; o, =1.0; o, =1.3. 
13 分 析 方 法 
空气 动力 学 粗糙 度 是 近 地 表 风速 为 零 的 高 度 ， 
是 反映 地 表 对 风阻 抗 的 重要 参数 ” ,粗糙 度 越 大 ， 
表明 地 面 对 风 的 削弱 能 力 越 强 。 其 计算 公式 如 下 : 
Z,=exp|(V,InZ,-V,InZ,/(V,-V,)] 0) 
式 中 : Z 为 空气 动力 学 粗糙 度 (em); Vi V, 分 别 是 
Z, Z, Wi es BEY UE ms) 
植株 防风 效能 计算 公式 如 下 : 


n..=(1- “= }x 100% (8) 


2 


v: 


式 中 : n, N (x, 习 处 的 防风 效率 ; x 为 距 植株 的 水 
平 距离 (m); z 为 距 地 面 高 度 (m); .为 植株 前 后 
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(x, 处 的 风速 值 (ms '); v, 为 裸 沙 地 在 z 高 度 处 - (a) 风 洞 试验 = 
的 初始 风速 值 (ms )。 T 018 _。 白 刺 
1.4 风 洞 试验 概述 
为 验证 模拟 结果 的 可 靠 性 ,本 文 利 用 已 有 风 洞 5 0.12 
试验 中 验证 相关 结论 。 该 试验 在 中 国 林业 科学 院 入 oos 
沙漠 林业 实验 中 心 进行 ,采用 直流 开口 吹 气 式 风  & oo 
洞 , 宽 、 高 均 为 2 m, 试 验 段 为 30 m。 风 洞 通过 变频 a 
器 调节 风速 ,本 文 设计 风速 为 6ms'。 由 于 本 文 主 0.00 


要 利用 风 洞 试验 测定 的 各 株 型 灌 从 阻 沙 效益 来 验 
证 相关 模拟 结果 ,因此 ,以 下 内 容 主 要 详 述 了 有 关 
风 洞 试验 中 测定 植株 阻 沙 效益 的 试验 布置 。 采 用 
阶梯 式 集 沙 仪 测定 灌 从 周围 不 同 高 度 层 的 输 沙 
量 。 选 择 具有 代表 性 的 白 刺 、 梭 梭 、 沙 拐 束 灌 从 ,在 
样 方 内 算 选 高 lm 左右 的 同龄 单 株 汐 木 ,并 将 具有 
代表 性 的 植株 取 回 试验 地 ,精确 测定 其 植被 特征 参 
数 。 将 采集 到 的 3 类 生长 季 灌 丛 单 株 固定 在 风 洞 试 
验 区 的 中 轴线 处 ,在 株 后 1 日 处 布设 一 个 1m 高 .50 
层 的 阶梯 式 积 沙 仪 , 每 层 入 口 截面 为 2 cmx2 cm, 
为 确保 沙 源 充足 ,在 试验 段 布设 5 cm FRAUD AR, Ub 
样 经 70 目 得 网 过 滤 。 每 种 汶 丛 单 株 在 6m's 风速 
下 进行 吹 蚀 试 验 , 吹 蚀 时 间 为 2 min。 吹 蚀 结束 后 取 
出 各 高 度 层 的 集 沙 槽 ,用 精度 为 干 分 之 一 的 天 平 测 
定 各 个 高 度 层 的 积 沙 量 ,并 计算 输 沙 率 。 每 类 灌 从 
试验 重复 3 次 ,并 在 每 次 吹 蚀 试验 结束 后 补充 沙 量 
并 平整 沙 面 。 同 时 设置 空白 对 照 组 , 即 在 试验 段 不 


放任 何 植株 的 条 件 下 , 测 得 输 沙 率 作为 对 照 ,对 比 
分 析 各 灌 从 的 阻 沙 效应 。 
2 结果 与 分 析 


2.1 风沙 两 相 流 场合 理性 验证 

为 验证 本 文风 沙 流 场 设置 的 合理 性 ,将 模拟 结 
果 与 已 有 风 洞 试验 结果 所 进行 对 比 。 研 究 表明 风沙 
运动 中 沙 粒 体积 分 数 与 输 沙 率 的 分 布 状况 一 致 。 
因此 ,本 文 将 风 洞 试验 中 6 ms'( 即 摩 阻 风速 为 
0.3752 ms ') 风 速 下 , 株 后 1H 处 白 刺 和 裸 沙 地 的 输 
沙 率 数据 作为 对 照 ,分 析 同 条 件 下 数值 模拟 所 得 白 
刺 沙 粒 体 积分 数 随 高 度 的 分 布 值 ,结果 如 图 2 所 
示 。 株 后 1 五 处 输 沙 率 与 沙 粒 体 积分 数 随 高 度 变化 
规律 大 体 一 致 , 均 表 现 为 随 高 度 层 的 增加 逐渐 减 
小 ,总 体 呈 现 对 数 减 小 趋势 ,此 规律 符合 沙 粒 运 移 
特征 和 风沙 流 结构 特点 ,说 明 模拟 结果 具有 和 较 高 的 
可 靠 性 ,可 以 较为 准确 地 反映 植株 周围 的 流 场 
特征 。 


0.05 0.10 0.15 0.20 0.25 0.30 0.35 
高 度 /m 

7 (b) 数值 模拟 

6 一 白 刺 


沙 粒 体积 分 数 (x10) 


0.05 0.10 0.15 0.20 0.25 0.30 0.35 
高 度 /m 


图 2 输 沙 率 和 沙 粒 体积 分 数 随 高 度 分 布 


Distribution of sand transport rate and sand volume 


Fig. 2 
fraction with height 


2.2 单 株 灌木 的 防风 效果 

2.2.1 植株 周围 风速 流 场 特征 比较 ”为 了 对 比 3 种 
株 型 植被 对 流 场 的 影响 ,选择 2m 高度 处 风速 为 6 
ms! 的 风速 廊 线 进行 模拟 (下 文 速度 均 与 此 相同 ， 
不 再 获 述 ), 图 3 展示 了 3 种 株 型 灌 从 的 风速 分 布 特 
征 。 由 图 3 可 知 ,在 相同 来 流 条 件 下 ,3 种 植株 上 方 
均 出 现 了 加 速 区 ,迎风 侧 出 现 遇 阻 减速 区 ,同时 埋 
风 侧 出 现 了 率 流 减速 区 、 涡 旋 加 速 区 以 及 恢复 区 。 
气流 在 株 前 -2 HH 处 出 现 扰动 ,形成 遇 阻 减速 区 ,其 
原因 为 灌 从 枝 型 较为 密集 ,枝条 从 生 增 加 受 风 面 
R ,加速 了 风能 的 损耗 。 气 流 到 达 植 株 前 方 时 ,一 
部 分 气流 汇聚 抬升 并 在 植株 上 方形 成 加 速 区 ; 另 一 
部 分 气流 穿 过 植株 ,经 过 涡 旋 作用 ,在 植株 背风 侧 
风速 急剧 下 降 , 从 而 形成 了 亲 流 减速 区 ,该 区 是 植 
株 的 主要 防护 区 ,其 大 小 直接 影响 植株 的 防护 效 
果 。 株 后 , 随 着 株距 的 增加 ,风速 逐渐 恢复 至 正常 
水 平 ,从 而 植株 失去 防护 能 力 。 值 得 注意 的 是 , 植 
株 周围 出 现 了 明显 的 涡 旋 加 速 区 ,3 类 植株 涡 旋 加 
速 区 的 范围 大 小 表现 为 沙 拐 更 > $ > 白 刺 ,而 3 
者 对 应 枝 下 高 度 依次 为 0.12 m .0.08 m,0.03 m, 表 明 
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图 3 灌木 周围 风速 云图 
Fig.3 Cloud chart of wind speed around shrubs 


涡 旋 加 速 区 呈现 随 枝 下 高 度 增 加 而 明显 增 大 的 趋 
势 ,与 已 有 研究 结论 相似 2 。 

2.2.2 株 后 涡流 分 析 已 有 人 研究 表明 ,植株 防 沙 的 
气动 性 原理 主要 是 局 部 风速 的 降低 和 气流 的 循环 
泡 ( 即 涡 旋 流 ) 7 。 图 4 展示 了 不 同 株 型 灌 丛 周围 的 
风速 矢量 。 从 图 中 可 以 看 出 背风 侧 主要 存在 3 个 涡 
流 , 其 中 涡流 [I 是 主要 防护 区 域 ,日 受 灌 从 株 型 影 
响 , 两 涡流 的 位 置 不 同 。 白 刺 、 梭 梭 及 沙 拐 束 的 涡 
流 工分 别 位 于 株 后 0.5~3.5 H、0.5~5 H、0.5~5.5 H 


处 ,涡流 开 的 位 置 依次 为 4~5 H、6~7 H.6.5~7.5 Ho 
分 析 3 个 涡流 的 形成 原因 如 下 :由 于 气流 穿 过 植株 
时 受到 枝叶 挤 压 ,越过 植株 后 过 流 断 面 突然 增 大 ， 
从 而 形成 向 上 向 下 运动 的 分 流 ,向 下 运动 的 部 分 
气流 在 距离 植株 不 远 处 抬升 ,又 返回 植株 方向 ,从 
而 形成 涡流 工 ,并 成 为 植被 积 沙 的 主要 区 域 ;涡流 开 
形态 虽 较 小 ,但 强度 较 弱 ,也 会 导致 大 量 沙 粒 在 此 
沉积 ,其 形成 与 涡流 I 的 下 沉 气 流 与 涡 旋 加 速 区 的 
气流 相 碰 撞 有 关 , 气 流 挤 压 后 形成 涡流 荆 ; 受 涡流 I、 
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图 4 灌木 周 目 
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Fig. 4 Vector diagram of wind speed around shrubs 


TAS Sg se , UTE FAR IE GI Tat I. ER BL 
沙 能 力 与 涡流 的 强度 及 尺度 有 直接 关系 ,涡流 范围 
越 大 强度 越 弱 , 沙 粒 越 易 沉积 ,从 而 植株 的 阻 沙 效 
益 越 好 。3 类 植株 涡流 I 范围 大 小 依次 为 梭 梭 > 沙 
HR > 白 刺 ORBEA > > FH. RR 
与 沙 拐 吏 受 涡 旋 加 速 区 的 影响 ,涡流 I 虽 范围 大 ,但 
强度 也 较 大 。 因 此 ,在 积 沙 初始 阶段 , 流 经 梭 梭 、 沙 
拐 囊 的 挟 沙 气流 主要 在 株 后 6~7 H 处 积 沙 , 而 白 刺 
主要 沉积 在 株 后 3 五 附近 。 

2.2.3 植株 后 方 风速 廊 线 对 比分 析 ”为 进一步 探究 
3 类 株 型 灌 从 在 垂直 方向 上 的 防护 效果 ,在 6 mes” 
风速 条 件 下 ,提取 3 类 株 型 株 后 1 互 处 垂直 截面 , 绘 
制 风 速 廊 线 图 ,并 以 裸 沙 地 为 对 照 。 根 据 图 5 可 知 ， 
3 类 株 型 的 风速 轮廓 线 存在 明显 差异 。 首 先 表现 在 
近 地 表 起 始 风 速 不 同 , 白 刺 的 起 始 风 速 最 小 ,为 
1.57 ms ', 梭 梭 与 沙 拐 束 的 起 始 风 速 均 大 于 裸 沙 
地 ,表明 在 株 后 1 处 的 近 地 表 区 域 , 白 刺 具有 抗 风 


2.5: r 


1 2 3 4 5 6 7 8 
风速 /(m*s-!) 


图 5 单 株 灌木 后 1 了 H 处 风速 廓 线 图 
Fig. 5 Wind speed profile of a single shrub 1 H later 


蚀 作用 ,而 梭 梭 `. 沙 拐 训 反而 会 促进 植株 周围 发 生 


风蚀 ,此 现象 值得 在 生态 恢复 工程 中 进一步 研究 。 
此 外 ,3 类 株 型 的 防护 高 度 也 存在 明显 差异 , 白 刺 、 梭 
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进一步 分 析 3 类 株 型 灌 从 与 裸 沙 地 的 风速 廊 线 
还 可 以 发 现 , 裸 沙 地 风速 轮廓 线 基本 呈现 对 数 形式 
增长 ,而 灌 丛 株 后 风速 轮廓 线 并 没有 形成 规律 的 对 
数 函 数 曲 线 , 而 是 随 高 度 增加 呈现 竖 向 “W” 变 化 的 
趋势 ,气流 垂 向 分 布 可 依据 其 变化 趋势 划 为 低空 减 
速 区 . 涡 旋 区 .高空 加 速 区 。 受 植株 的 阻挡 作用 , 白 
刺 风 速 廊 线 在 0.05~0.3 m 高 度 范围 内 逐渐 减 小 ,并 
在 0.3 m 高 度 处 达到 风速 极 小 值 ,风速 由 1.57 mss 
降低 至 0.18 mes, 风速 降幅 为 88.53% ,该 区 域 即 为 
低空 减速 区 ;随后 气流 在 0.3~0.8 m 高 度 层 内 进入 涡 
旋 区 ,风速 随 高 度 增加 呈现 先 小 幅度 增 大 后 减 小 的 
趋势 ;在 0.8 mm 高 度 处 以 后 ,风速 以 较 大 的 加 速度 呈 
现 快速 增长 的 趋势 ,植株 也 逐渐 失去 防护 能 力 ,此 
区 即 为 高 空 加 速 区 。 与 白 刺 不 同 , 梭 权 的 低空 减速 
区 出 现在 0.1~0.4m 高 度 层 内 ,风速 从 3.81 m- s~ KE 
低 至 0.18 ms ,风速 降幅 为 95.27%, 并 在 0.4 m 高 度 
出 现 风 速 极 小 值 ;随后 出 现 涡 旋 区 ,风速 在 0.4~0.8 
m 出 现 小 范围 波动 ,并 在 0.8 m 人 处 风速 达到 次 小 值 ; 
气流 在 0.8 m 以 上 呈现 随 高 度 增 加 而 增 大 , 即 高 空 
加 速 区 。 沙 拐 惠 具有 相似 的 变化 规律 ,其 低空 减速 
区 为 0.05~0.4 m, 风 速 降幅 为 64% ,0.4~0.8 m 为 涡 旋 
区 ,风速 出 现 先 增长 后 降低 的 特征 ,并 在 0.8 m 高 度 
出 现 风 速 极 小 值 ,随后 进入 高 空 加 速 区 。 总 体 而 
言 ,不 同 株 型 的 风速 轮廓 线 存在 差异 ,风速 极 小 值 
出 现 高 度 不 同 , 各 植株 最 优 防护 高 度 也 存在 明显 区 
别 , 白 刺 . 梭 梭 ` 沙 拐 束 的 最 优 防护 高 度 依次 为 0.2~ 
0.4 m,0.3~0.6 m,0.8~1 m。 值 得 注意 的 是 ,0.8 m 以 
上 高 度 ,各 植株 均 呈 现 风速 随 高 度 增加 而 增 大 的 趋 
势 , 表 明 0.8 m 人 处 是 风速 发 展 的 重要 分 界 点 ,此 现象 
值得 进一步 探究 。 
2.24 株 后 空气 动力 学 粗糙 度 分 析 由 公式 (7) 可 
知 ,观测 高 度 和 风速 的 选择 是 计算 空气 动力 学 粗糙 
度 的 关键 。 根 据 前 文 分 析 可 知 ,0.5~2 m 间 风速 轮廓 
线 的 变化 可 以 代表 整个 灌 从 层 的 风速 特征 , 故 选择 
3 类 株 型 植被 背风 侧 和 裸 沙 地 0.5 m 和 2 m 两 个 高 度 
的 风速 值 , 并 计算 得 出 3 类 株 型 株 后 不 同 株距 处 的 
空气 动力 学 粗糙 度 ,计算 结果 如 图 6 所 示 。 根 据 图 6 
可 知 , 随 株距 增 大 ,3 类 株 型 株 后 的 粗糙 度 逐 渐 减 
小 ,并 在 15 HH 处 与 裸 沙 地 粗 烟 度 基 本 持平 。3 类 株 
型 在 株 后 1 处 的 粗糙 度 最 大 , 白 刺 、 梭 梭 、 沙 拐 囊 
的 粗糙 度 依 次 是 裸 沙 地 的 26 倍 .20 倍 .13 倍 。 白 刺 
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图 6 灌 从 株 后 空气 动力 学 粗糙 度 


Fig.6 Aerodynamic roughness of shrub plants 


的 空气 动力 学 粗糙 度 明 显 高 于 其 他 两 种 株 型 , 沙 拐 
束 的 粗糙 度 最 小 ,表明 白 刺 削弱 近 地 表 风 力 的 能 
更 强 , 具 有 更 好 的 防风 作用 。 
2.2.5 植株 防风 效能 分 析 分别 提取 6 ms 风速 条 
件 下 各 高 度 层 不 同 株距 下 的 风速 值 ,通过 公式 (8) 
计算 单 株 治 木 的 防风 效能 。 分 析 图 7, 株 前 水 平 防 
风 效 益 随 株距 减 小 而 增 大 ,尤其 在 -2~0 五 范围 内 ， 
防风 效益 剧 增 ; 垂 向 防风 效益 随 高 度 增加 而 减 小 ， 
但 不 存在 负 值 ,表明 各 灌 从 在 其 高 度 范围 内 对 风速 
均 具 有 明显 的 前 弱 作 用 。 株 后 根据 各 高 度 层 防风 
效益 的 变化 趋势 ,可 将 其 分 为 近 地 表 区 (0.1 m、0.2 
m,0.3 m) 和 中 高 空 区 (0.5 m、0.7 m .0.9 m)。 在 近 地 
表 区 , 株 后 水 平 防风 效率 随 株距 增加 呈现 波浪 式 变 
化 ; 垂 向 防风 效益 在 株 后 5 也 范围 内 ,呈现 随 高 度 增 
加 而 增 大 的 趋势 。 在 中 高 空 区 的 垂 向 防风 效益 上 ， 
白 刺 、 梭 梭 防 风 效 益 随 高 度 增加 呈现 波动 式 降低 ， 
而 沙 拐 囊 呈 现 波动 式 增 大 ;在 水 平方 向 上 整体 呈现 
随 株距 增加 而 降低 的 趋势 , 白 刺 、 梭 梭 、 沙 拐 枣 的 防 
风 效 益 分 别 在 12 H、11 H,13 H 处 基本 下 降 至 50%， 
且 在 该 区 内 , 株 后 1H 处 防风 效果 最 好 。 以 沙 拐 束 
为 例 ,3 个 高 度 层 最 大 防风 效益 依次 为 89.8% 、 
98.1% .94.49% , 均 出 现在 株 后 1 HAE. 

进一步 分 析 图 7 可 知 , 除 梭 梭 、 沙 拐 束 在 0.1 m 
高 度 层 防风 效益 出 现 负 值 外 ,各 株 型 在 低 于 其 高 度 
内 均 产 生 了 有 效 的 防护 ,但 各 灌 从 的 最 优 防护 区 域 
存在 明显 差异 。 白 刺 最 优 防护 区 域 出 现在 0.3 m 高 
REJE PRIA 5 HH 范围 内 ,此 区 域内 植株 防风 效益 均 大 
于 90%; 梭 梭 最 优 防 护 区 域 出 现在 0.5 m 高 度 层 PR 
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图 7 灌木 植株 防风 效益 
Fig. 7 Protection benefit of shrub plants 


Ja 4 H WEN , in > ip UE LE 0.9 m 高 度 层 、 株 
后 2~5 H 范 围 内 。 总 体 而 言 ,各 株 型 在 -2~10 HD 
围 内 均 可 有 效 降低 风速 ,但 3 类 株 型 灌 从 的 区 域 防 
风 效 益 显 著 不 同 , 白 刺 在 近 地 表 区 具有 较 好 的 防风 
效果 ,而 梭 梭 、 沙 拐 束 在 植株 中 高 区 域 具有 较 好 的 


2.3 单 株 灌木 周围 沙 粒 运动 特征 分 析 
2.3.1 植株 前 后 输 沙特 征 分 析 ” 当 风 速达 到 起 沙 风 
速 时 ,地 表 沙 粒 就 会 在 风力 作用 下 开始 以 蠕 移 ER 
移 . 悬 移 的 形式 运动 ,进而 形成 风沙 流 。 灌 从 能 够 
拦 蕾 风沙 流 中 的 沙 粒 ,使 其 沉积 在 植株 周围 ,从 而 
形成 相应 的 灌 从 沙丘 ,但 不 同 株 型 灌 从 的 阻 沙 能 
不 同 。 为 了 探究 3 种 典型 株 型 荒漠 灌 从 的 阻 沙 作 
用 ,图 8 给 出 了 3 类 株 型 植株 在 6 ms"' 风 速 条 件 下 ， 
T=4 s、7=10 s 的 积 沙 云图 。 白 刺 前 期 迎风 坡 积 沙 量 
较 少 , 沙 粒 主要 沉积 在 背风 坡 , 但 与 植株 保持 一 定 
距离 ;随时 间 增 加 ,背风 坡 积 沙 点 逐渐 和 问 植 株 方向 
靠近 ,7=10s 时 , 株 后 积 沙 点 已 从 0.5 五 处 蔓延 至 植 
株 底部 , 且 背 风 坡 积 沙 长 度 由 3 再 变化 至 5 H, 株 前 
积 沙 长 度 由 0.5 于 增加 至 3.5 也 ,此 时 株 前 PR AFE 
沙 长 度 为 8.5 H, PRP URE AR BRU tt Be > HE 
要 集中 在 植株 附近 ,背风 坡 积 沙 量 较 多 , 积 沙 距离 
同样 表现 出 随时 间 增 加 逐渐 靠近 植株 的 趋势 ,7=10 
s 时, 株 后 积 沙 点 已 从 2 理 处 到 达 1H 处 , 株 后 积 沙 
长 度 大 致 从 4 H 变 化 至 5 HR KAY BEM cm 
增长 至 6 cm,10s 时 总 积 沙 长 度 为 6H。 沙 拐 杰 迎风 
坡 积 沙 最 少 ,背风 坡 在 4s 时 存在 较 大 的 积 沙 距离 ， 
表现 为 在 距 植 株 2.5 H 处 积 沙 , 随 时间 增 加 至 10 s, 
株 后 积 沙 点 到 达 1 HH 处 ,总 积 沙 长 度 增加 量 为 1.5 
H, 积 沙 高 度 增加 量 为 5 cm,10s 时 总 积 沙 长 度 为 
4.5 日, 相 比 其 他 株 型 , 沙 拐 束 的 阻 沙 效益 相对 较 差 。 
2.3.2 植株 前 后 地 表 蚀 积 特征 分 析 在 6m:s 的 风 
速 条 件 下 ,气流 与 沙 床 相 互 作用 ,撞击 和 挤 压 使 得 
沙 床上 的 沙 粒 以 多 种 形式 运动 ,从 而 导致 沙 床上 的 
沙 粒 发 生 重 排列 ,进而 使 沙 床 堆积 率 发 生变 化 。 在 
实际 工作 中 , 沙 床 堆积 率 不 易 精 准 检 测 ,但 数值 模 
拟 可 以 解决 这 一 难点 。 为 了 探究 不 同形 态 灌 从 植 
株 前 后 的 地 表 蚀 积 状况 ,本文 初始 条 件 设 置 了 厚度 
0.05 mm 上 且 均匀 的 沙 床 ,单一 粒 径 初始 堆积 率 为 
0.625 ,通过 提取 7T=10s 时 刻 的 治 程 沙 床 堆积 率 , 反 
映 不 同 株 型 对 区 域 蚀 积 状况 的 影响 。 根 据 植 株 周 
围 风沙 运 移 规 律 以 及 前 人 相关 人 研究 ,文中 提出 假 
设 : 若 0.05 m 以 下 高 度 层 沙 粒 体积 分 数 低 于 0.625， 
且 同 时 满足 0.05 m 以 上 高 度 层 沙 粒 体 积分 数 在 0.02 
附近 , 则 表明 该 区 域 发 生 风 蚀 ; 若 0.05 m 以 上 高 度 
层 体积 分 数 均 大 于 0.02, 且 此 刻 0.05 em 以 下 高 度 层 
的 沙 粒 体积 分 数 在 0.625 附近 , 则 说 明 该 区 域 出 现 
风 积 现象 。 

从 图 9a 可 以 明显 看 出 ,风沙 流 在 距 入 口 3 m Ze 
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8 灌 从 周转 积 沙 云图 


Fig. 8 Sand accumulation cloud atlas around shrubs 


右 处 受到 白 刺 影响 ,0.05 m 及 以 下 高 度 层 沙 粒 体积 
分 数 始终 保持 在 0.625 左右 ,0.05 em 以 上 高 度 层 沙 
粒 堆积 率 随 株 距 缩 短 而 增加 ,表明 株 前 发 生 了 风 积 
作用 。 在 白 刺 周围 ,0.05 cm 以 上 高 度 层 沙 床 堆积 率 
发 生 轻微 扰动 ,但 其 堆积 率 始终 高 于 0.5 ,表明 植株 
附近 发 生 了 风 积 作用 。 随 后 气流 逐渐 恢复 ,但 受 涡 
流 区 工 \、 卫 的 影响 形成 弱 风 人 蚀 区 、 风 积 区 。 梭 梭 和 沙 
拐 玉 在 株 前 一 定 范围 内 也 存在 不 同 程 度 的 风 积 作 
用 ,但 与 白 刺 不 同 的 是 ,两 者 在 植株 下 方 或 株 后 一 
定 范围 内 ,高 空 沙 粒 堆积 率 虽 增加 ,但 梭 梭 在 株 下 
0.04 m 高 度 层 WHETER F 0.04 m、0.05 m 高度 
层 , 沙 粒 体积 分 数 均 出 现 低 于 0.02 的 状况 ,表明 该 
区 域 发 生 风 蚀 作 用 ,风蚀 区 域 范 围 分 别 距 人口 5~ 
5.5 m、4.5~6 m; 气 流 穿越 植株 后 ,两 者 分 别 在 株 后 
距 入 口 约 6~10 m、6~9 m 发 生 风 积 作 用 。 总 体 而 言 ， 
3 类 株 型 灌 丛 对 地 表 蚀 积 微 地 貌 的 影响 存在 显著 差 
异 , 株 前 、 株 后 方 均 存 在 一 定 程度 的 风 积 状 况 ,但 除 
白 刺 外 ,其 他 两 种 植株 周围 出 现 不 同 程度 的 风蚀 
现象 。 


3 讨论 


3.1 不 同 株 型 灌 丛 防风 能 力 差异 的 原因 

沙 生 汐 丛 的 形态 差异 是 不 同 灌 丛 对 空间 IG 
水 等 资源 的 适应 策略 , 灌 丛 在 适应 不 同 生存 环境 时 
的 反应 特征 主要 体现 在 其 外 部 形态 上 "”。 生 存 环 
境 影响 植株 形态 的 发 育 进 程 ,而 植株 形态 也 会 影响 
其 自身 的 防风 阻 沙 效益 ”。 本 文通 过 模拟 发 现 , 单 
株 灌 从 周围 的 流 场 可 以 划分 为 5 个 区 域 :过 阻 减速 
区 、 上 方 加 速 区 、 闪 流 减 速 区 、 涡 旋 加 速 区 、 恢 复 
区 。 其 中 , 涡 旋 加 速 区 的 形成 与 植株 校 下 高 度 相 
关 。 较 大 的 枝 下 高 度 会 使 风沙 流 穿 过 时 ,附近 流 场 
产生 "* 狭 管 效 应 ””。 各 株 型 灌 丛 的 枝 下 高 度 不 同 ， 
导致 涡 旋 加 速 区 范围 存在 明显 差异 ,并 呈现 随 枝 下 
高 度 增 加 而 增 大 的 趋势 。 各 株 型 背风 侧 均 存 在 3 个 
涡流 ,涡流 强度 影响 沙 粒 的 沉积 位 置 。 由 于 坛 形 灌 
从 的 枝叶 密集 ,迎风 面 透风 系数 较 小 ,导致 涡流 I 
强度 较 小 , 即 在 风 影 区 风速 降幅 较 大 ,并 在 风 影 区 
形成 风 影 堆积 。 因 此 ,在 积 沙 初始 阶段 , 白 刺 主要 
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图 9 灌木 周围 蚀 积 状况 


Fig.9 Erosion around shrubs 


在 株 后 3 晴 处 积 沙 。 梭 形 、 帅 形 植株 底部 枝叶 分 布 
稀 玻 , 玻 透 度 较 大 ,气流 经 过 植株 孔 险 时 受到 挤 压 
而 加 速 ,形成 强度 较 大 的 涡流 工 。 相 比 涡 流 工 , 涡 
Dit E AY a BEBE, ATT SE BOI ie TE FR ie I) LE 
株 后 6~7 HH 处 积 沙 。 

3 类 灌 从 株 后 1 日 处 的 风速 廊 线 图 呈现 竖 问 


“W” 变 化 趋势 ,气流 垂 向 分 布 可 划分 为 低空 减速 区 、 
涡 旋 区 、 高 空 加 速 区 ,其 分 布 规律 与 植株 形态 相 
关 。 受 枝 下 高 度 的 影响 ,起 始 风 速 较 大 ,后 受 植株 、 
沙 粒 蚂 移 \ 跃 移 等 因素 影响 ,风速 随 高 度 增加 而 减 
小 ,从 而 出 现 低空 减速 区 ; 涡 旋 区 的 出 现 与 植株 本 
身 的 复杂 校 系 结 构 相 关 , 使 气流 穿 过 复杂 孔隙 时 出 
现 小 尺度 汕 流 ;高 空 加 速 区 则 是 受 植株 阻挡 影响 ， 
风速 随 高 度 增加 而 增 大 。3 类 株 型 灌 从 株 后 1 日 处 
的 风速 极 小 值 依次 出 现在 0.3 m、0.4 m、0.8 m 高 度 
处 ,最 优 防护 高 度 依次 为 0.2~0.4 m、0.3~0.6 m、0.8~ 
1 mm, 不同 株 型 的 防护 范围 及 最 优 防护 高 度 均 存在 
显著 差异 。 研 究 表明 ,防风 效应 取决 于 迎风 面 侧 影 
面积 的 大 小 ,植株 形态 决定 其 侧 影 面积 大 小 及 分 
布 ,从 而 影响 其 防护 范围 及 最 优 防护 高 度 ”。 白 
刺 形 态 低 矮 ,呈现 上 窜 下 宽 的 坛 形 , 近 地 表 区 域 枝 
系 繁 皮 ,最 大 侧 影 面积 出 现在 0.1~0.3 m 高 度 层 ， 
此 白 刺 可 以 有 效 降 低 近 地 表 区 风速 ,其 最 优 防护 高 
度 为 0.2~0.4 m。 梭 梭 整 体形 态 呈 现 中 间 宽 两 头 罕 
的 梭 形 ,具有 一 定 的 枝 下 高 度 , 植 株 底 部 形成 涡 旋 
加 速 区 ,从 而 容易 导致 局 部 风速 变 大 ,但 由 于 其 迎 
风 面 侧 影 面积 较 大 , 且 最 大 侧 影 面 积 出 现在 0.3~0.6 
m 高度 层 ,因此 ,其 具有 较 好 的 防风 效果 ,最 优 防风 
高 度 出 现在 植株 中 间 高 度 处 。 沙 拐 囊 呈现 上 宽 下 
窗 的 蝴 形 ,具有 较 大 的 枝 下 高 度 ,导致 其 在 近 地 表 
的 防风 能 力 较 弱 , 且 容易 引起 植株 根部 周围 的 侵 
刨 ,但 由 于 它 形 态 高 大 ,最 大 侧 影 面积 出 现在 0.6~1 
m 高 度 ,导致 其 在 0.6 m 以 上 高 度 层 具 有 较 好 的 防风 
效果 ,并 在 0.8~1 m 表 现 出 最 优 防护 高 度 。 

总 体 而 言 ,各 株 型 灌 从 的 防风 能 力 存 在 显著 差 
异 。 坛 形 植株 具有 较为 紧密 的 校 叶 结构 , 旦 其 底部 
侧 影 面积 较 大 ,可 以 有 效 降 低 近 地 表 风 速 ; 梭 形 ar 
形 植 株 虽然 底部 枝叶 相对 稀 玖 ,容易 形成 “ 狭 管 交 
应 ”, 但 两 者 形态 高 大 , 冠 层 中 高 部 分 透风 系数 小 ， 
因此 ,其 在 中 高 区 域 具有 较 好 的 防风 效益 。 
3.2 不 同 株 型 灌 从 阻 沙 能 力 差异 的 原因 

3 类 株 型 灌 丛 的 阻 沙 能 力 呈 现 白 刺 > RH > 沙 
PE , 总 积 沙 长 度 依次 为 8.5 H、6H、4.5 Ho 与 防风 
效益 相似 , 株 型 也 会 影响 植株 的 阻 沙 能 力 。 已 有 研 
究 表明 ,半球 形 冠 形 植株 比 锥 体 冠 形 及 梭 形 冠 形 植 
株 的 截 沙 能 力 更 强 吕 。 本 文通 过 研究 发 现 白 刺 的 
积 沙 范围 最 大 , 梭 梭 次 之 , 沙 毛囊 最小, 且 检 权 ` 沙 
拐 现在 植株 周 于 存在 明显 风蚀 区 ,因此 , 相 较 于 梭 
R 沙 拐 囊 , 白 刺 具有 更 好 的 防风 阻 沙 能 力 ,此 结论 
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与 已 有 试验 结果 相似 “””*。 其 原因 为 风沙 活动 主要 
发 生 在 近 地 表 0~0.3 mm 处 ,而 白 刺 呈现 下 宽 上 罕 的 
坛 形 ,贴近 地 面 生长 , 当 风 沙 流 经 白 刺 时 ,底部 枝 系 
结构 紧密 ,对 沙 粒 有 明显 的 阻挡 作用 ,导致 越过 和 
透 过 植株 的 沙 粒 减少 ,起 到 良好 的 阻 沙 作用 ,与 其 
容易 形成 沙 堆 的 特点 一 致 下 。 沉 积 在 3 类 灌 从 周围 
的 积 沙 量 始终 表现 出 背风 侧 高 于 迎风 侧 , 产 生 上 述 
积 沙 迁移 规律 的 主要 原因 为 :植株 前 后 对 风沙 流 的 
削减 能 力 存 在 差异 ,风沙 流 流 至 植株 时 ,迎风 侧 风 
速 昌 有 所 降低 ,但 大 部 分 气流 受 挤 压 抬 升 ,此 时 部 
分 动能 转化 为 重力 势能 ,加 之 风沙 流 流 经 植株 内 部 
孔隙 时 与 植株 相 撞 ,损失 的 动能 较 大 ,气流 挟 沙 能 
力 进一步 降低 ,因此 ,更 多 沙 粒 堆积 在 背风 坡 , 其 至 
形成 小 沙 堆 。 

就 灌 从 对 地 表 蚀 积 微 地 貌 影 响 的 研究 中 ,本 研 
究 发 现 3 类 株 型 灌 从 株 前 、 株 后 方 均 存在 一 定 程 度 
的 风 积 状况 ,但 梭 梭 、 沙 拐 束 分 别 在 距 入 口 5~5.5 m, 
4.5~6 m 处 存在 不 同 程度 的 风蚀 现象 ,这 与 两 者 存在 
较 大 的 枝 下 高 度 有 关 。 白 刺 株 前 部 分 风沙 流 从 植 
株 周 围 加 速 通 过 , 带 走 部 分 沙 粒 ,但 由 于 白 刺 近 地 
表 区 冠 层 侧 影 面积 大 ,导致 植株 下 方 沙 粒 沉 积 量 远 
大 于 吹 蚀 量 ,因此 ,其 株 下 并 未 发 生 风 蚀 而 是 出 现 
沙 粒 沉 积 。 而 流 经 权 梭 . 沙 拐 束 的 风沙 流 受 “ 狭 管 
效应 ”的 影响 ,在 植株 下 方 加 速 通过 并 带 走 大 部 分 
沙 粒 ,说明 两 类 灌 丛 根部 及 周围 存在 一 定 的 风蚀 
现象 。 

在 土壤 风蚀 防治 研究 中 ,研究 单一 灌木 或 灌 丛 
防护 效益 的 同时 ,还 要 考虑 灌木 林带 的 群体 防护 效 
应 和 机 制 兰 ?%。 合 理 配 置 灌 从 可 以 使 水 分 条 件 限 制 
下 的 有 限 植 被 达到 良好 的 防治 风蚀 效果 。 利 用 数 
值 模拟 研究 单一 植株 周围 流 场 的 目的 之 一 是 设计 
更 有 效 的 稀 玻 灌 从 防护林, 如 Guo 等 史 通 过 数值 模 
拟 发 现行 距 在 2.5~3 H 的 沙 苹 防 护林 防护 效益 最 
好 ,而 柳树 防风 林 的 最 优 间距 为 1~1.5 H。 本 研究 发 
现 白 刺 、 梭 梭 、 沙 拐 玉 的 最 优 防护 高 度 依次 为 0.2~ 
0.4 m,0.3~0.6 m、0.8~1 m, 而 3 者 的 阻 沙 能 力 依 次 为 
白 刺 > RIR > 沙 拐 枣 。 综 上 所 述 ,在 实际 风蚀 防治 
工程 中 ,可 采用 白 刺 与 梭 梭 、 沙 拐 玉 相 结合 的 配置 
方式 , 既 能 发 挥 白 刺 的 阻 沙 作用 ,又 能 利用 梭 梭 、 沙 
拐 东 在 中 高 空 区 较 大 的 防风 效应 。 


4 结论 


本 文 利用 欧 拉 双流 体 非 定常 模型 ,对 3 种 典型 
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株 型 荒漠 灌 从 周围 的 风沙 流 场 进行 了 人 研究 ,并 通过 
已 有 风 洞 试验 验证 了 流 场合 理性 ,得 出 以 下 结论 : 

(1) 单 株 治 从 周围 的 流 场 可 分 为 5 个 区 域 ,其 
中 , 涡 旋 加 速 区 呈现 随 枝 下 高 度 增加 而 增 大 的 趋 
势 。 株 后 均 出 现 3 个 涡流 ,但 同一 涡流 的 强度 存在 
差异 ,导致 在 积 沙 初始 阶段 , 梭 形 、 帅 形 植株 主要 在 
株 后 6~7 HH 处 积 沙 , 而 坛 形 主要 在 株 后 3 日 处 积 沙 。 

(2) 3 类 株 型 株 后 1 日 处 的 风速 轮廓 线 呈 现 竖 
向 “W” 变 化 趋势 。 受 植株 最 大 侧 影 面积 高 度 层 的 影 
响 , 坛 形 、 梭 形 \ 帅 形 的 风速 极 小 值 依次 出 现在 0.3 
m、0.4 m、0.8 m 高 度 处 ,最 优 防 护 高 度 依次 为 0.2~ 
0.4 m、0.3~0.6 m、0.8~1 m。3 类 株 型 株 后 的 空气 动 
力学 粗 烟 度 逐渐 减 小 , 且 坛 形 植株 的 粗糙 度 明显 高 
于 其 他 两 种 株 型 。 

(3) 3 类 株 型 在 -2~10 H 范 围 内 均 可 有 效 降 低 
风速 , 株 后 近 地 表 区 防风 效益 表现 为 坛 形 > 梭 形 > 
帅 形 ;中 高 空 区 防风 效益 均 随 株 距 增加 而 减 小 。3 
类 株 型 最 优 防护 范围 存在 明显 差异 , 坛 形 植株 最 优 
防护 区 域 出 现在 0.1~0.3 m 高 度 、 株 后 5 HH 范围 内 ; 
梭 形 最 优 防护 区 域 在 0.3~0.6 m 高 度 、 株 后 4H 范 围 
内 ;而 在 0.6~1 m、 株 后 2~5 HH 范围 内 蚌 形 防风 效益 
均 高 于 90%。 

(4) 在 7=10s 时 , 灌 从 株 前 、 株 后 方 均 存在 一 定 
程度 的 风 积 状况 , 白 刺 、 梭 梭 、 沙 拐 束 周围 总 积 沙 长 
度 分 别 为 8.5 H.6H\4.5 H。 受 “ 狭 管 效 应 ”的 影响 ， 
梭 梭 、 沙 拐 束 分 别 在 距 入 口 5~5.5 m 、4.5~6 m 处 存在 
不 同 程度 的 风蚀 现象 。 相 比 而 言 , 白 刺 的 阻 沙 效果 
较 好 ,但 在 实际 防风 固沙 工程 建设 中 应 使 其 与 梭 
梭 、 沙 拐 束 结合 , 既 能 发 挥 白 刺 的 阻 沙 作用 ,又 能 利 
用 梭 梭 、 沙 拐 束 在 中 高 区 较 大 的 防风 效果 。 
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Numerical simulation of the influence of typical shrub types on 
wind-sand flow field 


YAN Qing’, LI Juyan’, YIN Zhongdong', LIU Jinmiao', LIU Hongcai' 
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Abstract: The purpose of this study was to provide a theoretical basis for the rational selection of different 
vegetation types for wind and sand control in arid and semiarid areas. Fluent software was used to numerically 
simulate the flow field around three types of typical strains (altar-shaped, shuttle-shaped, and broom-shaped) of 
scrub to analyze the influence of different plant forms on wind and sand flow and verify the results using existing 
wind tunnel tests. Results showed that (1) The flow field around the three types of plants can be divided into five 
zones, and three eddies exist behind the plants. Due to the intensity of the eddies, during the initial phase of sand 
accumulation, shuttle- shaped and broom- shaped plants accumulated sand primarily at 6-7 H after the plant, 
whereas altar-shaped plants accumulated sand at 3 H. (2) Affected by the height layer of the maximum profile 
area of the plant, the minimum wind speed at 1 H after the three types of plant shrubs appeared at heights of 0.3, 
0.4 m, and 0.8 m, and the optimal protection range of height was 0.2- 0.4 m, 0.3- 0.6 m, and 0.8- 1 m, 
respectively. The aerodynamic roughness of the three plant types decreased gradually, and the roughness of the 
altar-shaped plants was significantly higher than that of the other two plant types. (3) All three types of strains 
could effectively reduce wind speed in the range of -2- 10 H. The wind protection benefits in the near-surface 
area after the strains are altar-shaped > shuttle-shaped > broom-shaped the wind protection benefits in the mid- 
altitude area all decrease with increasing plant distance. (4) At T= 10 s, the total duration of sand accumulation 
around the three plants was 8.5 H, 6 H, and 4.5 H, respectively, and wind erosion existed to different degrees at 5- 
5.5 m and 4.5-6 m from the entrance for Haloxylon ammodendron and Calligonum mongolicum, respectively. 
Compared with other plants, Nitraria tangutorum exerted a better sand-blocking effect. Hence, it is recommended 
to combine N. tangutorum with H. ammodendron and C. mongolicum in the construction of wind and sand 
fixation projects, so that the sand- blocking property of N. tangutorum can be effectuated, and the better wind- 
blocking effect of H. ammodendron and C. mongolicum can be utilized at medium and high altitudes. 
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